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In this article we show that the asymptotic iteration method (AIM) allows one to numerically 
find the quasinormal modes of Schwarzschild and Schwarzschild de Sitter (SdS) black holes. An 
added benefit of the method is that it can also be used to calculate the Schwarzschild anti-de Sitter 
(SAdS) quasinormal modes for the case of spin zero perturbations. We also discuss an improved 
version of the AIM, more suitable for numerical implementation. 
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I. INTRODUCTION 



The study of quasinormal modes (QNMs) of the Schwarzschild black hole is an old and well established subject 
(for a recent review see reference [l| and references therein), where the various frequencies have been well determined, 
typically by applying a Frobenius series solution approach, leading to continued fractions for QNM boundary con- 
ditions, a la Leaver [2j. Recently a new method for obtaining analytic/numerical solutions of second order ordinary 
differential equations with bound potentials has been developed called the asymptotic iteration method (AIM) 0], 
which was found to be closely connected to continued fractions developed from exact WKB solutions, see reference 
^ , [i[ and those listed in therein. 1 

In this article we will demonstrate that the AIM can also be applied to the case of black hole QNMs, which have 
unbounded (scattering) like potentials. On a related note, the AIM was used to find QNMs for Scarf II (upside-down 
■ Poschl- Teller- like) potentials Q , based on observations made by one of the current authors @ relating QNMs from 
quasi-exactly solvable models. Indeed bound state Poschl- Teller potentials have been used for QNM approximations 
previously by inverting black hole potentials 0- However, the AIM docs not require any inversion of the black hole 
' potential as we shall show. 

In this paper we shall focus on spherically symmetric backgrounds in four dimensions with field equations of the 
!* ■ form: 

^ + - v(x)] = o , (i) 



where V(x) is a master potential of the form [l[ 



V(r) = f(r) 

and dx = dr/f(r) where 



, n 2s(2M (4~s 2 )A 
1 1 V r 3 6 



(2) 



, . 2M A , 

f(r) = 1 - — - 3 r 2 (3) 
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1 The phrase "continued fraction" used for WKB solutions should not be confused with the type of continued fractions developed from 
Frobenius series for QNMs by Leaver 0. In the rest of the paper "continued fraction" will mean those generated from Frobenius series. 
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with cosmological constant A. Here s = 0,1,2 denotes the spin of the perturbation: scalar, electromagnetic and 
gravitational, (for half-integer spin see references In the following we shall use the definition that QNMs are 

defined as solutions of the above equations with boundary conditions: 

"> ( e— z ^ -co , ^ 

for e~ lut time dependence (which corresponds to ingoing waves at the horizon and outgoing waves at infinity). We 
shall see how the AIM can not only be applied to Schwarzschild and Schwarzschild de Sitter (SdS) QNMs for general 
spin, s, but also for the Schwarzschild anti-de Sitter case (SAdS), at least for spin zero (s = 0) fields. 

The structure of the paper is as follows. In the next section we describe how the AIM can be applied to exterior 
eigenvalue problems for the Schwarzschild background. Then in Section ITTll we show how the AIM can be applied to 
backgrounds with more than one horizon, such as SdS. In Section ITVl we show how the AIM can be applied to scalar 
fields on SAdS, while in Section [V] we discuss the results and make conclusions. In an appendix we present a brief 
derivation of the continued fraction method for SdS for spin s = 1, 2 fields. 



II. SCHWARZSCHILD BLACK HOLES 



To explain the asymptotic iteration method (AIM) we shall start with the simplest case of the radial component 
of a perturbation of the Schwarzschild metric outside the event horizon . For an asymptotically flat Schwarzschild 
solution (A — 0): 



where from dx = dr/f(r) we have 



/M-l-^, (5) 
r 



x(r)=r + 2Mln(^-l) , (6) 



for the tortoise coordinate, x. Note, for the Schwarzschild background the maximum for this potential, in terms of r, 
is given by 

ro = ^j^(t(t + 1) - - + + I) 2 + + 1)- + - 2 ) 1/2 ) , (7) 

where a = 1 — s 2 [ll| ■ 

The choice of coordinates is somewhat arbitrary and in the next section (for SdS) we will see how an alternative 
choice leads to a simpler solution. First, consider the change of variable: 

2M 

£ = 1 -— > ( g ) 
with < £ < 1. In terms of £, equation ([l} then becomes, 

d 2 ip l-3£ dip 



de c(i - o da 



4MV /(/ + !) 1-* 2 1^ = . (9) 



£ 2 (i-0 4 £(i-£) 2 £(i-0 



To accommodate the out-going wave boundary condition ip — > e lLJX = e lu >( r + 2Min ( r / 2M !)) as (x,r) — > oo in terms 
of £ (which is the limit £ — > 1) and the regular singularity at the event horizon (£ — > 0), we define 

m = r aj(« (1 _ e)- 2 '^ e W x g) , (10) 

where the Coulomb power law is included in the asymptotic behaviour (cf. reference 0] equation (5)). The radial 
equation then takes the form: 

X" = A (0x' + *o(0x , (11) 

where 

Ao(0 = €(l-0 2 ' 

16Af 2 o; 2 (£ - 2) - 8Miu(l - Q +1(1+ 1) + (1 - s 2 )(l -Q M „. 
*o(0 = g(l-fl 2 ' ( } 
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Note that primes of \ denote derivatives with respect to £. 

The crucial observation in the AIM is that differentiating the above equation n times with respect to £, leaves a 
symmetric form for the right hand side: 

X {n+2) = K(OX + Sn(Z)X , (14) 

where 

K(0 = + sn-i(0 + and «n(0 = *Ui(0 + *o(0*»-i(0 ■ (is) 

For sufficiently large n the asymptotic aspect of the "method" is introduced, that is: 

s n(£) s n-l(£) _ af j-\ 

Al = Vl =f) ' (16) 
where the quasinormal modes are obtained from the "quantization condition" 

S n = s n X n -i - s„-iX n = , (f 7) 

which is equivalent to imposing a termination to the number of iterations (l2j . This leads to the general solution 



x(0 = cx p 



(18) 



for integration constants C\ and C2, which can be determined by an appropriate choice of normalisation. Note, that 
for the generation of exact solutions C\ = 0. 



Improved AIM 

One unappealing feature of the recursion relations in equations (|I5[) is that at each iteration one must take the 
derivative of the s and A of the previous iteration. This can slow down the numerical implementation of the AIM and 
also lead to problems with numerical precision. To circumvent these issues we have developed an improved version 
of the AIM which bypasses the need to take derivatives at each step. This greatly improves both the accuracy and 
speed of the method. We expand the A„ and s n in a Taylor series around the point which the AIM is performed, £: 

00 

An(0 - XXCs-O', (19) 

2=0 
OO 

«n(f) = XX^-O*, (2°) 
i=0 

where the d n and are the i-th Taylor coefficient's of A„(£) and s n (£) respectively. Substituting these expressions 
into equations [T51 leads a set of recursion relations for the coefficients: 

i 

< = (* + lK^+^- 1 +5>oC i i J (21) 

k=0 

i 

di = (i + lRt 1 1 +^4C- fe i- (22) 

In terms of these coefficients the "quantization condition" (fTT|) can be reexpressed: 

- d° n _ lC l = (23) 

and thus we have reduced the AIM into a set of recursion relations which no longer requires the derivative operator. 

Observing that the right hand side of equations (|21[) and (j2"2"j) involve terms of order at most n — 1, one can recurse 
these equations until only c\ and d l terms remain (that is the coefficients of Aq and sq only). However for large 
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TABLE I: Quasinormal modes to 4 decimal places for gravitational perturbations (a — —3) where the fifth column is taken 
from reference fO /. Note that the imaginary part of the n — 0, £ = 2 result in JlH] has been corrected to agree with reference 
0/ (see [*]). Note, if the number of iteration in the AIM is increased to say about 50 then we find agreement with reference 0/ 
accurate to 6 significant figures. 



£ 


n 


^Leaver 


cjaim (after 15 iterations) 


U)WKB 


2 





0.3737 - 0.0896 i [*] 


0.3737 - 0.0896 i 
(<0.01%)(<0.01%) 


0.3732 - 0.0892 i 
(-0.13%) (0.44%) [*] 




1 


0.3467 - 0.2739 i 


0.3467 - 0.2739 i 
(<0.01%)(<0.01%) 


0.3460 - 0.2749 i 
(-0.20%) (-0.36%) 




2 


0.3011 - 0.4783 i 


0.3012 - 0.4785 i 
(0.03%) (-0.04%) 


0.3029 - 0.4711 i 
(0.60%)(1.5%) 




3 


0.2515 - 0.7051 i 


0.2523 - 0.7023 i 
(0.32%) (0.40%) 


0.2475 - 0.6703 i 
(-1.6%) (4.6%) 


3 





0.5994 - 0.0927 i 


0.5994 - 0.0927 i 
(<0.01%)(<0.01%) 


0.5993 - 0.0927 i 
(-0.02%) (0.0%) 




1 


0.5826 - 0.2813 i 


0.5826 - 0.2813 i 
(<0.01%)(<0.01%) 


0.5824 - 0.2814 i 
(-0.03%) (-0.04%) 




2 


0.5517 - 0.4791 i 


0.5517 - 0.4791 i 
(<0.01%)(<0.01%) 


0.5532 - 0.4767 i 
(0.27%) (0.50%) 




3 


0.5120 - 0.6903 i 


0.5120 - 0.6905 i 
(<0.01%)(-0.03%) 


0.5157 - 0.6774 i 
(0.72%)(1.9%) 




4 


0.4702 - 0.9156 i 


0.4715 - 0.9156 i 
(0.28%)(<0.01%) 


0.4711 - 0.8815 i 
(0.19%) (3.7%) 




5 


0.4314 - 1.152 i 


0.4360 - 1.147 i 
(1.07%) (0.43%) 


0.4189 - 1.088 i 
(-2.9%)(5.6%) 


4 





0.8092 - 0.0942 i 


0.8092 - 0.0942 i 
(<0.01%)(<0.01%) 


0.8091 - 0.0942 i 
(-0.01%) (0.0%) 




1 


0.7966 - 0.2843 i 


0.7966 - 0.2843 i 
(<0.01%)(<0.01%) 


0.7965 - 0.2844 i 
(-0.01%) (-0.04%) 




2 


0.7727 - 0.4799 i 


0.7727 - 0.4799 i 
(<0.01%)(<0.01%) 


0.7736 - 0.4790 i 
(0.12%) (0.19%) 




3 


0.7398 - 0.6839 i 


0.7398 - 0.6839 i 
(<0.01%)(<0.01%) 


0.7433 - 0.6783 i 
(0.47%) (0.82%) 




4 


0.7015 - 0.8982 i 


0.7014 - 0.8985 i 
(-0.01%) (-0.03%) 


0.7072 - 0.8813 i 
(0.81%)(1.9%) 



numbers of iterations, due to the large number of terms, such expressions become impractical to compute. We avert 
this combinatorial problem by beginning at the n — stage and calculating the n + 1 coefficients sequentially until 
the desired number of recursions is reached. Since the quantisation condition only requires the i = term, at each 
iteration n we only need to determine coefficients with i < N — n, where N is the maximum number of iterations to 
be performed. 

The QNMs that we calculate in this paper will be determined using this improved asymptotic iteration method. 
Indeed in a previous work of ours for spheroidal harmonics (l3j , we found that the AIM was somewhat slower than 
the continued fraction method. However, the Mathematica code implementing this new version of the AIM is found 
to be on a par with the continued fraction approach; these can be downloaded from (l4j . 



III. SCHWARZSCHILD dS 



The QNMs for Schwarzchild gravitational perturbations are presented in TableU] however, to further justify the use 
of this method, it is instructive to consider some more general cases. As such, we shall now consider the Schwarzschild 
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de Sitter (SdS) case, where we have the same WKB-like wave equation and potential as in equation ([T]), but now 

,.2 



2A7 , r" 
/(r) = l- — -A T , 



(24) 



where A > is the cosmological constant. Interestingly, the choice of coordinates we use here leads to a simpler 
AIM solution, because there is no Coulomb power law tail; however, in the limit A = we recover the Schwarzschild 
results. Note that although it is possible to find an expression for the maximum of the potential in equation ([1}, for 
the SdS case, it is the solution of a cubic equation, which for brevity we refrain from presenting here. In our AIM 
code flij we use a numerical routine to find the root to make the code more general. 

In the SdS case it is more convenient to change coordinates to £ = 1/r |15| . which leads to the following master 
equation (cf. equation 



d 2 ip p' dtp 



t{l + !) + (!- s 2 ) (2M£_ - (4 - s 2 ) A 



p 



p 



4' = o , 



where we have defined 

p = f - 2Alf - A/3 => 
It may be worth mentioning that for SdS we can express (l5| : 

e wx = (£ - £i)^(£ - £ 2 )^(£ - £ 3 )^ 



p' = 2£(1 - 3M0 



(25) 



(26) 



(27) 



in terms of the roots of f(r), where £i is the event horizon and £2 is the cosmological horizon (and K n is the surface 
gravity at each £„). This is useful for choosing the appropriate scaling behaviour for QNM boundary conditions. 

Based on the above equation an appropriate choice for QNMs is to scale out the divergent behaviour at the 
cosmological horizon: 2 



(28) 



which implies 



pu" + (p' - 2iu)u' - 



A 



l) + (l- s *)[2M£,-(A-s 2 ) — 



in terms of £. Furthermore, based on the scaling in equation 
implies: 



u = , (29) 
the correct QNM condition at the horizon £1 



u{x) = (£-£1) ^x(x) , 
where k± is the surface gravity at the event horizon £1: 



2dr 



= Mi{ - - 



1 A 



36 ' 



(30) 



(31) 



with £1 = 1/ri, where T\ is the smallest real solution of f(r) = 0, implying p — 0. The differential equation then 
takes the standard AIM form: 



X" = Ao(£)x' + so(£)x 



where 



Ao(£) 
*o(f) 



2ico 



«i(£-£i) 



2iu> 



2A/e-(4-, 2 )A 



■ioj 



«i(£-£i) 



flU) 



[p' - 2iu) 



«i(£-£i) 



(32) 

(33) 
. (34) 



2 Note that this is opposite to the case presented in reference [Ta |. where they define the QNMs as solutions with boundary conditions: 
ip(x) oc e^" x as x — ¥ ±00, cf. equation for e %LJt time dependence. 
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where in Table [TTI results are presented for SdS with s = 2. 

For completeness in the Appendix we derive a set of three-term recurrence relations for the continued fraction 
method, valid for electromagnetic and gravitational perturbations (s = 1,2). It may be worth mentioning that via the 
AIM we can treat the s = 0,1,2 perturbations on an equal footing (for the scalar case the continued fraction method 
reduces to a five-term recurrence relation, see the Appendix). 

TABLE II: Quasinormal modes to 6 significant figures for Schwarzschild de Sitter gravitational perturbations (s — 2) for 1 = 2 
and £ = 3 modes. We only present results for the AIM method, because the results are identical to those of the continued 
fraction method after a given number of iterations (in this case 50 iterations for both methods). The n — 1,2 modes can be 
compared with the results in (52/ for s — 2. 



A {£ 


= 2) 


n=\ 


n = 2 


n — 3 







0.373672 


- 0.0889623 i 


0.346711 


- 0.273915 i 


0.301050 


- 0.478281 i 




0.02 


0.338391 


- 0.0817564 i 


0.318759 


- 0.249197 i 


0.282732 


- 0.429484 i 




0.04 


0.298895 


- 0.0732967 i 


0.285841 


- 0.221724 i 


0.259992 


- 0.377092 i 




0.06 


0.253289 


- 0.0630425 i 


0.245742 


- 0.189791 i 


0.230076 


- 0.319157 i 




0.08 


0.197482 


- 0.0498773 i 


0.194115 


- 0.149787 i 


0.187120 


- 0.250257 i 




0.09 


0.162610 


- 0.0413665 i 


0.160789 


- 0.124152 i 


0.157042 


- 0.207117 i 




0.10 


0.117916 


- 0.0302105 i 


0.117243 - 


0.0906409 i 


0.115876 


- 0.151102 i 




0.11 


0.0372699 - 


0.00961565 i 


0.0372493 - 


0.0288470 i 


0.0372081 - 


0.0480784 i 


A [t 


= 3) 


n= 1 


n = 2 


n = 3 







0.599443 


- 0.0927030 i 


0.582644 


- 0.281298 i 


0.551685 


- 0.479093 i 




0.02 


0.543115 


- 0.0844957 i 


0.530744 


- 0.255363 i 


0.507015 


- 0.432059 i 




0.04 


0.480058 


- 0.0751464 i 


0.471658 


- 0.226395 i 


0.455011 


- 0.380773 i 




0.06 


0.407175 


- 0.0641396 i 


0.402171 


- 0.192807 i 


0.392053 


- 0.322769 i 




0.08 


0.317805 


- 0.0503821 i 


0.315495 


- 0.151249 i 


0.310803 


- 0.252450 i 




0.09 


0.261843 


- 0.0416439 i 


0.260572 


- 0.124969 i 


0.257998 


- 0.208412 i 




0.10 


0.189994 


- 0.0303145 i 


0.189517 - 


0.0909507 i 


0.188555 


- 0.151609 i 




0.11 


0.0600915 - 


0.00961888 i 


0.0600766 - 


0.0288567 i 


0.0600469 - 


0.0480945 i 



IV. SCHWARZSCHILD AdS 



There are various approaches to finding QNMs for the SAdS case (an eloquent discussion is given in the appendix 
of reference [HI). One approach is that of Horowitz & Hubeny [T^], which uses a series solution chosen to satisfy the 
SAdS QNM boundary conditions. This method can easily be applied to all perturbations (s = 0,1,2). The other 
approach is to use the Frobenius method of Leaver Q , but instead of developing a continued fraction, the series must 
satisfy a boundary condition, such as Dirichlct, at infinity fl5j . 

The AIM does not seem easy to apply to metrics where there is an asymptotically anti-de Sitter background, 
because for general spin, s, the potential at infinity is a constant and hence would include a combination of ingoing 
and outgoing waves, leading to a sinusoidal dependence [ijj]- However, for the scalar spin zero (s = 0) case, the 
potential actually blows up at infinity and is effectively a bound state problem. In this case the AIM can easily be 
applied as we show below. 

Let us consider the scalar wave equation in SAdS spacetime, where A = — 3/i? 2 , and R is the AdS radius. The 
master equation takes the same form as for the graviational case [1] except that the potential becomes: 

H-^)(H= 2(r " 1>(r3+ ; +2)(r3+1) - <»» 

Here for simplicity we have taken the AdS radius R = 1, the mass of the black hole M — 1, and the angular momentum 
number / = 0. Hence, the horizon radius r+ = 1. Thus, with this choice we can compare with the data in Table 3.2 
on page 37 of reference [l9[ (see Table |Hl] below) . 

To implement the AIM we first look at the asymptotic behavior of ip. One could obtain this by using the WKB 
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approximation: 



|^-vf /4 

As r — > r + = 1, the potential V goes to zero. In addition, 

rdr 



i 



(r-l)(r 2 +r + 2) 4 
4 , 

dr 



-ln(r - 1) + 



-ln(r - 1) 



±iu/4 



_4(r- 1) 

^ „ e ±<[* In(r-l)] „ (r _ 1)± «,/4 ~ f 1 - 

For QNMs we choose the "out-going" (into the black hole) boundary condition. That is, 



ip ~ e 



1 



—iui/i 



(36) 

(37) 
(38) 

(39) 

(40) 

(41) 



On the other extreme of our space, r — ¥ oo, the potential goes to infinity. This is a crucial difference from the case 
of gravitational perturbations. In that case, the potential goes to a constant. However, in the scalar case, as r — > oo, 



vy 



dr 



Vuj 2 - V 



1> 



1 

- H 

1\ 1/2 
r) + 

iV2 In 
r 



and to implement the Dirichlct boundary condition, we take 



For the AIM one possible choice of variables is 



1 



r 

and we see that to accommodate the asymptotic behaviour of the wavcfunction we should take 

Finally, after some calculation we find the scalar perturbation equation to be: 

x" = A (0x' + s (Ox , 

where 

-iuq + 2[-4 + 2(9 + 4>/2)f - (21 + 10^2)^ + 4(2 + v^ 3 ] 



A 



So 



2£g 

[ 4iw[9 + 8V2- 2(7 + 5%/2)£ + (6 + 4,V2)f]q + uj 2 (-1 + 2 (-40 + 41£ - 20£ 2 + 4£ 3 ) 
-4[4 - 5£ + 2^ 2 ] [-8(3 + 2\/2) + 8(10 + 7V2)£ - (91 + 64%/2)£ 2 + (34 + 24\/2)£ 3 



(42) 
(43) 
(44) 

(45) 

(46) 

(47) 

(48) 
(49) 

(50) 



(51) 



8 



and q = (—4 + 9£ — 7£ 2 + 2£ 3 ). Using the AIM we find the results presented in Table [TTT1 below, which are discussed 
in the next section. 



TABLE III: Comparison of the first few quasinormal modes to 6 significant figures for Schwarzschild anti de Sitter scalar 
perturbations (s = 0) for £ = modes with r+ = 1. The second column corresponds to data \19l] using the Horowitz & Hubeny 
method flU /. while the third column is for the AIM using 70 iterations. Note the mismatch for real part of the n — 3 mode in 
Jldlj . see [*] below; we have confirmed this using the Mathematica notebook provided in 0j/ . 



n 


HH method 


AIM 





2.7982 - 


2.6712 i 


2.79823 


-2.67121 i 


1 


4.75849 - 


5.03757 i 


4.75850 


-5.03757 i 


2 


6.71927 - 


7.39449 i 


6.71931 


-7.39450 i 


3 


8.68223[*] 


- 9.74852 i 


8.68233 


-9.74854 i 


4 


10.6467 - 


12.1012 i 


10.6469 


-12.1013 i 


5 


12.6121 - 


14.4533 i 


12.6125 


-14.4533 i 


6 


14.5782 - 


16.8049 i 


14.5788 


-16.8050 i 


7 


16.5449 - 


19.1562 i 


16.5457 


-19.1563 i 



V. RESULTS AND DISCUSSION 



We first applied the AIM to the Schwarzschild background defined in Section HH setting M = 1 and choosing £ = £o 
(the black hole maximum), which typically leads to the fastest converging results. In Table Q] the QNMs for graviton 
perturbations using the AIM are compared to some other methods. These results show that after fifteen iterations 
the AIM is in good agreement with Leaver's method to four significant figures, with the disagreement becoming most 
pronounced for the lowest £ modes with larger values of n. We also found that by increasing the number of iterations 
the AIM can be made to agree exactly with the results of Leaver's continued fraction method [H to any number of 
decimal places. 

For completeness we have also included results from an approximate semi-analytic 3rd order WKB method 
(more accurate semi-analytic results with better agreement to Leaver's method can be obtained by extending the 
WKB method to 6th order f20l|). In this regard, it may be worth mentioning that a different semi- analytic approach 



has recently been discussed by Dolan and Ottewill |21| , which has the added benefit of easily being extended to any 
order in their pcrturbative scheme. 

We also discussed spin s = 0, 1,2 in SdS (see Table ITT|) where results to six significant figures are presented for the 
graviton (s = 2) case. Identical results were generated by the AIM and continued fraction method, both after 50 
iterations. Though results are presented for n = 1, 2 and 3; £ = 2 and 3 modes only, the AIM is robust enough to be 
applied to any other case; where like for the A = case, agreement with other methods in more extreme parameter 
choices would only require further iterations. For comparison with the AIM we have also used the continued fraction 
method, and for completeness we presented a derivation in the appendix (for spin s = 1,2). 

As far as we are aware, only reference [22| (who used a semi-analytic WKB a ppr oach) has presented tables for 
general spin fields for the SdS case. We have also compared our results to those in [22( for the s = 0, 1 cases and find 
identical results (to a given accuracy in the WKB method). As briefly discussed in the appendix another benefit of 
the AIM is that no Gaussian Elimination procedure is required; unlike for continued fractions on Reissner-Nordtrom 
[23j or higher-dimensional Schwarzschild backgrounds [24|, for example. All that is necessary is to factor out the 
correct asymptotic behaviour at the horizon(s) and infinity (we showed this for higher-dimensional angular spheroids 
in [Hj]). Of course, there is no need to find the maximum of the black hole potential for continued fractions, unlike 
with the AIM. 

The reader might wonder about approximate results for cases where Poschl- Teller approximations can be used, such 
as for Schwarzschild and SdS backgrounds, e.g., see @, HH. I n fact when the black hole potential can be modeled 
by a Scarf like potential the AIM can be used to find the eigenvalues exactly [B[ and hence the approximate QNMs. 
Indeed the authors of reference || also independently suggested the idea of using the AIM to find black hole QNMs 
from exact potentials numerically (based on previous comments by one of the present authors Q). On a related note 
the large £ eikonial limit (e-g^ for A = see [ll|) is also easily confirmed numerically as we have verified using the 
continued fraction method [14| . 

The SAdS case is a little more subtle, and we briefly discussed some preliminary results for a spin zero field (the 
I = 0, r + = 1 case) in Table IIIII In the scalar case the WKB form of the potential is essentially a bound state 
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potential (it blows up at infinity) and allowed us to apply the AIM numerically. For other spins the potential goes 
to a constant at infinity and thus would be a combination of ingoing and outgoing waves. This implies a sinusoidal 
boundary condition at infinity (e.g., see [l8[) which is a little awkward to use with the AIM and somewhat out of the 
scope of this article. 

In summary, we have demonstrated how the AIM can also be applied to radial QNMs and not just to spheroidal 
eigenvalue problems [HI, [l3| . Given the fact that the AIM can be used in both the radial and angular wave equations 
[1 31 ] we expect no problems in obtaining QNMs for Kerr or Kerr-dS black holes for example. Furthermore, a perusal 
of the work on the AIM (e.g., see references [1, H, [l2[ and references therein) indicates that there are still some 
mathematical problems/questions of numerical refinement (such as the improved AIM we discussed), which was 
another reason for bringing this method to light. 

As for future work, it remains to be seen if the AIM can be tailored to handle asymptotic QNMs (see reference [25[ 
for an adapted version of the continued fraction method). However, given the close relationship between the AIM and 
the exact WKB approach [J it might be possible to adapt the AIM to find asymptotic QNMs [2(| [27} numerically or 
even semi-analytically. Also, there are still issues for the SAdS case, but we hope to report on more general results 
for this in the near future. 
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Appendix: Continued fraction method for SdS with s = 1,2 

For completeness and for the purpose of making an accurate comparison with the AIM method, for the SdS case 
we present some the steps that lead to the continued fraction method. As far as we are aware the actual recursion 
coefficients themselves do not seem to appear in the literature (e.g., see [HI). 

The correct scaling for QNM boundary conditions in SdS implies 

i>(x) = e i » x {Z-Z 1 )-%x(x) , (A.l) 

where we can expand x as a power series: 

X = J]a n y n , y = ^r => £ = (6-£i)i/ + 6, (A.2) 

where ao = 1. Note that the form of y, above, depends on (£2 — £1) i n order to give a converging solution in the 
Frobenius solution. It then follows that if we define ip = X^^Li a nJ/" ~ then 

1>(x) = e iu *{Z 2 -Z 1 )-%<p{x), (A.3) 

and hence the (£2 — £1) term factors out of the master equation ([23]). The rescaled equation for <p then takes an 
identical form to equation (|29|) : 



;2Af£ 3 -(4-s 2 )^ 



<p = , (A.4) 



where for general integer spin we multiplied both sides of the equation by £ 2 . We could then expand ip as a power 
series where 



•r 



' = £<*»«-- y"-* "'(ft - 6)- 1 , */ = $>n-^ n-^-l W-^- 2 (6 - a)" 2 , (A.5) 



n—1 x * n— 1 



where we should also power expand ^ 2 p and ^ 2 p' in a power series about £1 (for a detailed discussion of this point see 
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It is straightforward to see that this will lead to a five-term recurrence relation, because 
p(£i) = 0, for the u" term. This would then require the use of two Gaussian Eliminations [2 



2 p ~ 0(y 5 ), with 
| to reduce it to a 

three-term recurrence relation. Note, in the Reissner-N6rdtrom-(A)dS case, even the axial graviton perturbations 
reduce to a four-term recurrence relation [l6j . 

Thus, in order to obtain a set of three-term recurrence relations (which are simpler to deal with as no Gaussian 
Eliminations are required) we consider only spin s = 1,2. This drops the last constant term in equation (|A.4[) and 
allows us to divide out by a factor of £ 2 . Then, because p ~ 0(y 3 ) for the u" term, this leads to a three-term 
recurrence relation. After substituting these expansions into equation (|A.4[) and equating coefficients we obtain a set 
of recurrence coefficients: 



2(6 - Zi)-\n - - + 1) \{n - - + 1)6(1 - 3M6) 

K 



-£(£ + 1) - 2M(1 - s 2 )6 + (1 - 6Affi)(n )(n + 1) 

-2M(Z 2 -Z 1 )(n--+s)(n---s) , 



(A.6) 



which satisfy the following recurrence relations: 

a ai + P a Q = 
a n a n+ i + f3 n a n + 7„a„_i = . 

The exterior eigenvalue problem (QNMs) then leads to a continued fraction 

«o7i "172 



= /?o 



(A.7) 
(A.8) 

(A.9) 



where this equation can be solved for numerically and allows us to check the results with the AIM, see [14] . In all cases 
we find perfect agreement between the two results (to a given precision) for roughly the same number of iterations in 
both the AIM and continued fraction. Note that in the recurrence coefficients above, the overall factors of (£i — £2) 
cancel out in the continued fraction. 
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